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We describe an asymptotic procedure for deriving continuum equations from the kinetic theory 
of a simple gas. As in the works of Hilbert, of Chapman and of Enskog, we expand in the mean 
flight time of the constituent particles of the gas, but we do not adopt the Chapman-Enskog device 
of simplifying the formulae at each order by using results from previous orders. In this way, we are 

O 

able to derive a new set of fluid dynamical equations from kinetic theory, as we illustrate here for 
the relaxation model for monatomic gases. We obtain a stress tensor that contains a dynamical 
pressure term (or bulk viscosity) that is process-dependent and our heat current depends on the 

X. 

gradients of both temperature and density. On account of these features, the equations apply to a 
greater range of Knudsen number (the ratio of mean free path to macroscopic scale) than do the 
Navier-Stokes equations, as we see in the accompanying paper. In the limit of vanishing Knudsen 
number, our equations reduce to the usual Navier-Stokes equations with no bulk viscosity. 

PACS numbers: 05.20 Dd, 47.45 -n, 51.10 +y, 51.20 +d 
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I. INTRODUCTION 



The derivation of fluid equations from kinetic theory is often effected by an expansion 
of the solution of the kinetic equation in the Knudsen number, the ratio of the mean free 
path to the characteristic macroscopic scale. The method generally in use was developed 
by Chapman and Enskog II] whose aim in part was to remedy deficiencies in the earlier 



work of Hilbert El ID|. Chapman and Enskog derived the Euler equations in zeroth order 



and the Navier-Stokes equations at first order, starting from the Boltzmann equation of 
kinetic theory. But the Navier-Stokes equations do not adequately describe the dynamics of 
fluids when the Knudsen number is not very small. The inadequacies of the N-S equations 
for dealing with the problems of rarefied media are well documented and we shall cite 
empirical evidence for this when we compare our theory with experiment in papers II and 
III of this series. Problems arise for example in the study of shock waves since the thickness 
of a shock wave is generally of the order of the mean free path of the particles that make 
up the medium through which the shock propagates (paper III). However, the difficulty has 
little if anything to do with the nonlinearity of the shock waves, as we may see from the 
similar failure of the N-S equations to predict accurately the propagation of linear sound 
waves when their periods are comparable to the mean flight times of the constituent particles 
(paper II). 

To improve matters, one may try solving the kinetic equations directly [B3], but this 
leads to a problem in higher dimensions than is encountered in solving the continuum fluid 
equations. Alternatively, one may seek to derive from the kinetic theory a set of fluid 
equations with a greater domain of validity than the N-S equations. This was the aim of 
attempts to go to higher order in the development of the kinetic equation in Knudsen number. 
Higher order approximations, such as the Burnett or super Burnett developments || [| (and 
that of Woods f2"6fl ) have hardly improved matters and the resulting equations do not work 
well when the Knudsen number is not infinitesimal pT0[| . Moreover, the continuum equations 
found in these higher approximations are very complicated so that they do not seem to repay 
the effort involved in their use. 



Some hope for improvement of the situation was raised when Grad JTT| introduced his 
moment method, whose leading order results are the Navier-Stokes equations. However, the 



further development of Grad's method [13, IS] does not produce very rapid convergence nor 
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does it easily give very accurate solutions of the problems mentioned here (as we shall see in 



paper II). Levermore JT6| has proposed an alternate moment method that avoids the failures 
of causality that sometimes arise with the Grad method. The newer method has produced 
interesting results on shock structure [ff7[] but as yet has not given results with heat flux. 

Other approaches that are being explored are the resummation of the terms of the 
Chapman-Enskog expansion [T2I], [H], [Sj and the flux-limited diffusion theory jnj|; both 
of these approaches have proved useful in radiative transfer theory (the former in || |24| 
and the latter in [|l^]). A promising and interesting approach is by way of thermodynamic 
(or phenomenological) models [13|, |19| that give good agreement with experiment jn|. The 
relativistic extensions of these latter procedures ]12) produce hyperbolic systems and do not 
violate causality. 

Here we propose a modification of the asymptotic procedure of Chapman and Enskog. 
We begin, as they do, with an expansion in mean free paths of the constituent particles, a la 
Hilbert. However, we do not simplify the results at a given order by introducing expressions 
(mainly for time derivatives of fluid quantities) from previous orders as did Chapman and 
Enskog. As we shall see, this seemingly small difference in the methods makes for large 
differences in the results when they are applied to cases with finite mean free paths. As 
mentioned in [<J and more fully documented in paper II, even our first-order asymptotics 
provides acceptable accuracy when comparison is made with experimental results. 

In this paper, we derive continuum equations from the relaxation model of kinetic the- 
ory Jl|, [|, |25|] (also called the BGK model). In the relaxation model, the approach toward 
equilibrium takes place in a relaxation time, r, which is determined by the mean flight time 
of the constituent particles of the medium. So the analogue of Hilbert 's approach is here 
an expansion in relaxation time. In our version of this development, we shall not recycle 
the lower orders through the current order as is done in Chapment-Enskog theory. The 
fluid equations folowing from our procedure, generalize the Navier-Stokes equations, which 
we shall recover from ours by a simple development for the stress tensor and the heat flux 
in terms of r, which is effectively the Knudsen number. In that sense, a byproduct of our 
development is a relatively simple derivation of the Navier-Stokes equations, one that avoids 
some of the complications of the Chapman-Enskog method. 
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II. SOME KINETIC THEORY 



Consider a gas made up of identical particles of mass m obeying classical Hamiltonian 
dynamics. The phase space of a single particle is six-dimensional and its coordinates are 
spatial position, x, and velocity, v. The expected number of particles in a phase volume 
G?x<iv is /(x, v, t)dxdv where / is the density in phase space (or one-particle distribution 
function). For a Hamiltonian system the phase flow is incompressible in the sense that its 
six dimensional velocity is solenoidal, and the relaxation model of kinetic theory |T|, n, |2"5|| 
that we adopt here is 

V f = fojzl^ (2.1) 

T 

with the following notation. The streaming operator is 

V = d t + v%i + a% (2.2) 

where i = 1,2,3, repeated indices are summed, d x i = d/dx\ d v % = d/dv\ a 1 represents the 
acceleration of a particle caused by an external force and r is the time scale on which the 
system relaxes to the equilibrium, / - 

While we prefer to think of fl2.ip as a model, it may be also considered as an approximation 
to the Boltzmann equation where r and /o are taken as approximations to functionals of / 
that appear in the Boltzmann collision operator. With either interpretation, the relaxation 
time, r, is of the order of a collision time. Though r might in principle depend on v, we 
here assume that it depends on just the local values of the macroscopic fields. To give an 
expression for r we then need to specify the macroscopic quantities. 

The mass density of the fluid is defined as 

p(x,t) = f mfdv, (2.3) 



where the integration is over all of velocity space. We may also introduce the particle number 
density n = p/m. A second important macroscopic quantity is the mean, or drift, velocity, 
defined as 

u(x,t) = - J mvfdv . (2.4) 

We may then introduce the peculiar velocity c = v — u and use it to define the temperature 

as 



r=— /cVrfv, (2.5) 
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where R = k/m is the gas constant and is the Boltzmann constant. We assume throughout 
this work that / goes rapidly to zero as |v| — > oo, so that /-weighted integrals over velocity 
space are finite and well defined. 

The macroscopic quantities are all functions of x and t, and, when we speak of a temper- 
ature, we use the notion of local thermodynamic equilibrium, in which equilibrium formulae 
are used to describe nonequilibrium conditions locally in space and time. In the relaxation 
version of the kinetic equation, the interaction term on the right of ( |2.1|) , representing col- 
lisions amongst the particles, drives the system toward the equilibrium given by /q, which 
we take to be the local Maxwell-Boltzmann distribution in the frame locally comoving with 
the matter, namely 



This distribution depends on position and time only through its dependence on the local 
macroscopic quantities, n, u and T and it is moreover a local equilibrium solution of the 
Boltzmann equation. 

To make the description work well, we try to arrange that the fields on which /o depends 
are most nearly those of the real flow. This 'osculating' property of the assumed equilibrium 
may be imposed by matching conditions that are inherited from the theory of the Boltzmann 
equation. There the quantities that are conserved in two body collisions are the sums of the 
masses, the velocities and the energies of the of the two colliding particles. We let ipa with 
a = 0, 1, ...,4 be these summations invariants; that is — m , ^% — and ^4 = |mv 2 
with % = 1,2, 3. In Boltzmann theory, this conservation property results in the orthogonality 
of the ip a with the collision term in the kinetic theory. This property is generally ascribed 
to kinetic models including the relaxation model. Therefore when we multiply the right side 
of equation by any of the ip a and integrate over velocity space, we get zero. This gives us 
the matching conditions 



provided that we may assume that r does not depend on the velocity. This condition ensures 



Finally, we may specify that r should be of the order of the mean flight time of the 
constituent particles, that is, the mean free path of particles divided by the mean speed. 
If the collision cross-section of the particles is independent of velocity, then simple collision 






(2.7) 



that the true fluid fields are the ones appearing in (|2.6| ). 
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theory gives the result that r oc -^=, so that, generally speaking, r is generally a function 
of the macroscopic variables. 



III. EQUATIONS OF FLUID DYNAMICS 



When when we multiply (|2.1|) by rip a and integrate over v, the right side is zero because 
of the matching condition. If we add the assumption that r does not depend on v we obtain 
these macroscopic equations: 

d t p + V-{pu) = (3.1) 
p {d t u + u ■ Vu) + V ■ P = pa (3.2) 
^pR (d t T + u • VT) + P:Vu + V- Q = 0, (3.3) 

where the semicolon stands for a double dot product and we have introduced the pressure 
tensor (or second moment) 

P = J mcc/dv (3.4) 
and the heat flux vector (a third moment) given by 

Q = J lmc 2 c/rfv . (3.5) 

In deriving these equations, we have performed integrations by parts under the integrations 
over velocity space. These give no boundary terms if, as we assume, / goes rapidly enough 
to zero as |c| gets large. (Strictly speaking, the P in fl3.3|) should be the transpose, but as 
the stress tensor is symmetric in this work, we need not indicate this.) 

These macroscopic equations are a formal consequence of the kinetic equation, (|2.1| ). 
For them to be useful, we must supply adequate expressions for P and Q. The Navier- 
Stokes forms of these higher moments, which were first derived phenomenologically, were 
obtained for small mean free paths by the Chapman-Enskog method from the Boltzmann 
equation. The N-S equations have also been extracted from (|2 . 1| ) by this method ||. Those 
developments follow on the work of Hilbert, who introduced a series expansion in mean free 
path for / into the Boltzmann equation [|| |^, |j"0| . In the case of the relaxation model, 
matters are simpler because it has an essentially linear form and this permits a clearer 
statement of the underlying approach to deriving the fluid equations. The procedure is to 
develop in r, on the presumption that it is small compared to the macroscopic times that 



arise. Such an asymptotic development for the case of the relaxation model begins with an 
expansion of the form 

oo 

/ = E A™) 7 " • (3-6) 

m=0 

When we substitute ( |3.6| ) into ( |2.1| ) we obtain a series of equations for the /( TO ) . The first 
of these is the anticipated condition /( ) = fo- The second approximation is simply 

/(i) = -Vf • (3.7) 
When we introduce expression (|2.6| ) into this result we obtain the more explicit relation 



/(i) — ~fo 



( c 2 3\ 1 

V In p + )V In T H c ■ £>c 

H \2RT 2 J RT 



(3.8) 



Matters are quite simple up to this point but when we go to higher orders, we encounter 
terms arising from the derivatives of r. So we stop the expansion at this point since we 
already have enough in first order to obtain an interesting generalization of the equations of 
fluid mechanics. 

To use the asymptotic results for / in fl3.4|) and ( |3.5|) , we note that the higher moments 
can be written as 

P = P (0) + rP (1) + ... (3.9) 

and 

Q = Q(o) +rQ ( i) + ... (3.10) 

where 

P (n) = J mccf(n)dv and Q (n) = / imc 2 c/ ( „)dv . (3.11) 
Then we readily see that 

mc A f dc J eedtt (3.12) 

where c = |c| and e = c/c. 

With a little rearranging, the first of the two integrals in (|3.12 ) becomes proportional to 
a gamma function. The second can be written as 

eedVl=—I (3.13) 

where I is the unit tensor; its components are those of the Kronecker delta. Hence, we find 
that 

P (0) =pl (3.14) 



where 

p = RpT, (3.15) 

which we identify with the gas pressure. By contrast, the integral for the heat flux is odd in 
c and we find that Q(o) = 0. In short, the macroscopic equations at leading order are the 
Euler equations. 

To evaluate P(i), we note that 

D 



T> = — + c ■ V + a ■ d c 



where 



D 
Dt 



a + u • v 



and we use V and <9 X interchangeably. Then we have 

c 2 D 



/(i) - —fo 
where 



A + c • B + — — InT - (cc : V\nT)/(RT) + c • VlnT 

2RT Dt K Jn 1 2RT 



Dt-W B = Vln^ + (a-^)/(i ? T). 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



With the c-related factors thus in evidence, it is a straightforward matter to carry out 
the necessary integrals to evaluate Pm and Q(i). For this, we need the formula 



J e* e j e k e e dfl = ^| (5 ij 8 M + 5 ik 5 je + 5 ie 5 jk ) , 



(3.20) 



which may be verified by explicit evaluation of its components. If we include P(o), we may 
write the result of the integration as 



P 



DlnT 2„ 

• 3 V u 



I-fiE + 0{T 2 



where \i = rp, and 

For the heat current, we get 

Q = 



ox a oxj 3 



r/VT - r/TV lnp - ^ ( ^ - a ) I r 2 1 



(3.21) 



(3.22) 



(3.23) 



where r] = 

A more detailed derivation is given in |J. The expressions for P and Q involve not only 
the fluid dynamical fields, or lower moments of /, but their substantial derivatives as well. 
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Those derivatives are given by the fluid dynamical equations that, in turn, involve P and 
Q. So the higher moments are given only implicitly by ( |3.22| ) and ( |3.23| ) and we need to 
solve the fluid equations together with those equations to obtain explicit expressions for the 
higher moments. 



IV. DISCUSSION OF THE EQUATIONS 
A. The Entrance of Entropy 

Since the equations for P and Q are intertwined with the field equations themselves, we 
first rewrite these equations so as to clarify their meaning. We may introduce the continuity 
equation into ( |3.22| ) and so replace V • u by — p/p where the dot stands for D/Dt. The 
combination of T/T and of p/p that then appears suggests the introduction of the quantity 

For an ideal gas, this is the formula for the specific entropy, with C v = 3R/2 and 7 = 5/3. 
If we introduce S together with the definition of p, we may rewrite (|3.21|) as 

1 - l£. I _ + C(r 2 ) . (4.2) 



P 



Similarly, the momentum equation (|3.2|) can be used to rewrite the formula for the heat 



flux. Since ( |3.23|) contains Du/Dt, we may use ( |3.2|) to rewrite it as 



Q = - V VT - h -vV ■ I ^ I + rpE ] + G(r 2 ) , (4.3) 



'V 



where v = p/p. 

We may also use to convert ( |3.3|) into an evolution equation for the specific entropy. 
If we introduce ( |4.2j ) into that equation, we find a term pV • u so that p/p comes in by way of 
the continuity equation. The term in (|3.3| ) involving T/T combines with this and we obtain 
an equation for S. But S also appears in P and hence, when we gather the two apparitions 
of S together, we see that ( |3.3|) becomes 



2 DS 

pT(l - -rV • u)— = -pE : Vu - V • Q . (4.4) 

Though our equations may have an unfamiliar look, this is not because we have done 
anything unusual. Rather, we have omitted to do some things that are normally considered 
usual. So let us see how to get back to more familiar ground. 



B. The Euler and Navier-Stokes Equations 



We have obtained 



\i = rp and r\ = -rpR (4.5) 



and so we see that when r — > 0, we have P = pi + 0(r) and Q = 0(t), with S = 0{r) 
according to ([4.4|) . Thus we obtain the Euler equations in leading order together with 
entropy conservation. 



We see also that the term rS/C v in both (|4.2| ) and (|4.3Q is of order r 2 . Therefore, we 



find that 



P = pl- fiE + 0(r 2 ) 



(4.6) 



and 



Q = -r?VT + C(r 2 ) 



(4.7) 



when r is very small. Therefore, when the extra terms in our pressure tensor and heat flux 
are developed in r, we see that our forms differ from the standard Navier-Stokes terms in 
terms of order r 2 , which is perfectly allowable in first-order theories. 

In contrast to the conventional closure approximations, our expressions for P and Q 
depend on both the fluid fields and their derivatives. This means that these expressions 
must be solved in concert with the dynamical equations. To express P and Q explicitly in 
terms of the fluid fields, as in the usual closure relations, we would need to make expansions 
in r. As we have just seen, in the first order, we recover the Euler equations and in the 
second order we get the Navier-Stokes equations. Continuation of this development produces 
terms of all orders in r. Therein lies the crucial difference of our results from those of the 
Chapman-Enskog procedures. The extra terms in our development of the present results do 
not correspond to the higher theories based on Chapman-Enskog procedures that lead to 
the Burnett equations, as shall be explained in another place. All the terms in the present 
approximation come about from a first-order theory and the differences from standard theory 
arise in terms of second and higher order. When r is not infinitesimal, which it never is in 
practice, these terms do have an effect on the predictions of the theory. 
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C. Dynamical Pressure 



To clarify the meaning of the difference between our equations for a simple gas and those 
obtained with the Chapman-Enskog procedure we note that, as in |3.6| , we are writing the 
solution of the kinetic equation at any order as 

N 

f = f N + 0(T N+1 ) With f N =Y.hm)T m . (4.8) 

m=0 

From this we then obtain an approximation for the stress tensor in the form 

P = Pat + Mjv (4.9) 

where R N = 0(r N+1 ) is the error incurred in the truncation of the series. If we take the 
trace of Q4.9Q , we find 

tr P = tr Piv + tr Mjv • (4.10) 
As we see from the definition ( |3.4| ) of P, in the exact case, 

trP = 3p. (4.11) 

This result says that the total pressure has no contribution of a dynamical kind for a struc- 
tureless gas, as Maxwell and Boltzmann both recognized. The Chapman-Enskog procedure 
imposes this condition at every finite order so that 

tr¥% E = 3p. (4.12) 



When we impose (|4.12|) onto ( [4.10|) , we are forcing the requirement that tr Mat = 0, which 



overly constrains the results obtained by Chapman and Enskog and destroys any hope of 
improving convergence by mitigating their errors. However, all that we should demand of 
our successive approximations is that in the iVth approximation, 

ti¥ N = 3 P + 0{t n+1 ) (4.13) 

On introducing this less restrictive condition into ( [4.10[ ) we find that, in our procedure, 



tr R N = 0{t n+1 ) , (4.14) 

which does not incur the loss of generality that forcing the trace of Mat to vanish does. 
By keeping terms 0(r N+1 ) in our iVth approximation for P and Q, we leave open the 
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possibility of compensating for the errors caused by the truncation of the series for / by 
retaining suitable process-dependent effects of higher order. 

Thus we have in our present approximation a dynamical pressure in our approximate 
pressure tensor, F 1 = P( ) + tP(i), that is, 

= 3p + 0(r 2 ). (4.15) 

Our asymptotics suggests that the extra term compensates for the effect of truncation of 
the series for /. Since this term has contributions from all orders in r, it can in principle be 
very effective in extending the domain of validity of the theory. 

V. CONCLUSION 

We have illustrated our derivation of the fluid equations from kinetic equations by carrying 
out the procedure for the relaxation model of kinetic theory. The same procedure can be used 
on other forms of the kinetic equation. In the case of the Boltzmann equation, the procedure 
is quite similar, though the inversion of the linearized Boltzmann collision operator involves 
some technical issues that we shall take up elsewhere. However, by working out the case of 
the relaxation model we can more readily see the differences between our approach and the 
Chapman-Enskog procedure. 

In our derivation, we do not introduce slow times as in the Chapman-Enskog method. 
This means that we are not driven to expand the fluid variables (or slow quantities) in r 
as in C-E theory. Such expansions cause ambiguity in the application of initial conditions 
since it is not clear how to distribute the initial values over the various orders. Moreover, 
those expansions lead to a different sequence of approximations than ours. Chapman-Enskog 
theory gives the pressure tensor and the heat flux explicitly in terms of the fluid fields, a 
feature which results from invoking solvability conditions at each order. Our results do not 
produce explicit formulae for P and Q in terms of the fluid fields; rather, these quantities are 
expressed in terms of the fluid fields and their derivatives. Those derivatives appear in the 
field equations themselves so that we do produce a closed system of equations. Moreover, 
since the expansion variable r appears in the equations, we may further expand the equations 
to develop explicit formulas for P and Q accurate to any prescribed order. 

As we saw, the leading terms in the development of the formula for P and Q give us 
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tr Pi = 3p 



successively the Euler and Navier-Stokes equations, but the development need not stop 
there. That is, our finite formulae implicitly contain terms of all orders in r and this means 
that we may hope that they will produce high accuracy even when r is not infinitesimal. In 
any event, it is clear that we may expect a divergence between results from our system and 
those from the Navier-Stokes equations when r is not very small. 

J.B. Keller (private communication) has remarked that two theories with the same nomi- 
nal accuracy may have different domains of validity. In forcing the trace of Pa? to be exactly 
3p, order by order, the Chapman-Enskog method renounces the extra generality allowed by 
the freedom to choose higher order terms in an advantageous way. For us, the problem has 
been to select the best way to allow for the higher order corrections when trying to extend 
the domain of validity of the theory. We shall use as our test of validity the comparison 
with experiment given in the two following papers of this series. 
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